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Summary 


This paper describes the formulation of an elastic/viscoplastic constitu- 
tive model which was used to predict the measured behavior of graphite/ 
thermoplastic and graphitc/bismaleimide composite materials at elevated 
temperature. The model incorporates the concepts of overstress and effec- 
tive stress/strain to provide a simple formulation which was able to account 
for material behavior under monotonic tension or compression loads over a 
temperature range of 23° C to 200°C. Observed behavior such as stress re- 
laxation and steady state creep, in off-axis tension and compression tests, 
were predicted by the model. Material constants required by the model were 
extracted from simple off-axis test data. 


NOMENCLATURE 


A - quasistatic clastic/plastic material constant 
a 66 ■ potential function material constant 
E - elastic Young’s modulus 

F - time dependent inelastic strain rate function 
f - potential function 
G - elastic shear modulus 
II - overstress 

h(0) - function defined in equation 20 

K - elastic/viscoplastic material constant 

m - elastic/viscoplastic material constant 

n - quasistatic elastic/plastic material constant 

Q - stiffness matrix 

S - compliance matrix 

e - strain 

i - strain rate 

c - effective strain 

e - effective strain rate 

a - stress 

a - stress rate 

a - effective stress 

a - effective stress rate 

<7* - quasistatic stress 

a* - effective quasistatic stress 

0 - fiber angle relative to loading axis 

A - proportionality constant 

$ - overstress function 

Subscripts 

ij,k,l = 1 - (i, lamina principal directions 
x - laminate axial loading direction 
Superscripts 
e - elastic 

qp - quasistatic plastic 
p - plastic 
in - inelastic 
vp - viscoplastic 
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INTRODUCTION 


An elastic /viscoplastic constitutive model for polymer matrix composites was 
developed in order to meet two researcli objectives. The first objective was 
to provide analytical and experimental procedures for comparing the relative 
rate dependent response of several composite material systems under condi- 
tions of elevated temperature and high loads. The second objective was to 
help develop the analytical foundation for predicting various rate dependent 
phenomenon such as creep, stress relaxation and strain rate sensitivity in a 
laminated structure. 

In order to address these objectives, an elastic/viscoplastic constitutive 
model which incorporates the concepts of overstress and effective stress/strain 
was formulated. The model assumed that the lower bound of the rate de- 
pendent stress/strain behavior is described by rate independent (quasistatic) 
expressions. The rate dependent strain was decomposed into two terms to 
account for both loading and unloading. Material constants required by the 
model were extracted form simple off-axis tension or compression test data 
over a range of test temperatures. 


CONSTITUTIVE MODEL DESCRIPTION 

A variety of constitutive models have been suggested in recent years to model 
nonlinear, rate dependent stress/strain behavior in advanced composite ma- 
terial systems. A recent survey of several of these models and their ap- 
plications to polymer matrix composites can be found in Gates [1]. These 
models have been examined to determine their suitability for predicting the 
nonlinear, rate dependent stress/strain behavior of several advanced polymer 
matrix composites (PMC’s) at elevated temperature. In order to provide a 
means of analytically modeling the observed material behavior, the author 
has developer a modified form of the time dependent clastic/viscoplastic con- 
stitutive model given by Gates and Sun [2]. The quasistatic elastic/plastic 
constitutive model given by Chen and Sun [3] has also been used in the model 
development. This paper presents the mathematical formulation, assump- 


i 


3 


tions and possible shortcomings of the model. In addition, model predictions 
are compared to experimental results for two advanced PMC materials. 


Rate Dependent Elastic/ Viscoplastic Behavior 

For the elastic/viscoplastic behavior, the time dependent strain is assumed 
to be composed of elastic and viscoplastic components. For the general mul- 
tiaxial case, this is written as: 


(,j = 

r +c ,,p 

>j ' tj 


(i) 

where we have the individual constitutive relations: 



{U = \S) e {a} or 

% = 

elastic 

(2) 

and 

{*}-[SP{*} o'’ 

% = 

viscoplastic 

(3) 

If we further decompose the viscoplastic component into time dependent 
plastic and inelastic components, we can write': 


<i L ij i l ij 


and this implies the constitutive relations: 


=> 



and 



(4) 

( 5 ) 


It is noted that for this formulation, the superscript (p) refers to a rate 
dependent term while the rate independent, or quasistatic plastic, term is 
designated by the superscript (<//>). 

Using a procedure similar to that used in rate independent plasticity, a 
form of the associated flow rule for the rate dependent plastic strain can be 
written as: 



( 6 ) 


where A is a proportionality factor. Using the formulation given in Chen 
and Sun [3] and Kenaga [4], the potential function (/), which accounts for 
material anisotropy, is formed by assuming elastic behavior along the fiber 
direction and plane stress conditions. 


2 f{(Tij) = + 2 


(7) 
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Where cr 2 2 and a }2 arc the inplane transverse and shear stress components 
respectively. The single material constant is given by the a 66 term and can be 
found from test results on ofT-axis specimens. It is noted that as a GG decreases, 
a higher shear stress to transverse stress ratio is required to satisfy the yield 
condition. 

Using the overstress concept, the effective plastic strain rate is written as 
a function of overstress given by: 


F =<$(//)> ( 8 ) 

where H is the ” overs tress” and is defined as: 


II = (a - a*) (9) 

where (<r) is the rate dependent stress and (a*) is the rate independent or 
quasistatic stress. Using the potential function above, the effective stress is 
defined by: 

» = v/VK) (io) 

and a* is the effective quasistatic stress. A power law is assumed for the 
effective stress, effective quasistatic plastic strain relation. Functionally, this 
can be written as: 

i qp = A{a*) n (11) 

where A and n are material constants found from fitting the power law to the 
effective stress, effective plastic strain data. The Macaulay (<>) brackets 
imply a conditional statement which can be written in a general sense as 


< <&(//) >= 


<D (//) if // > 0 
0 if II < 0 


( 12 ) 


The concept of overstress and its relationship to viscoplastic strain in 
isotropic mctallics has been attributed to Malvern[5] and his work on high 
strain rate conditions during wave propagation. Additional references to 
overstress and its use in constructing viscoplastic models can he found in the 
work of Eisenberg and Yen[6,7] and Krempl and IIong[8]. 

The effective stress quantity allows for the dependence of the state of 
stress on the angle between the load and fiber directions and can be repre- 
sented in a general sense by equation (10). Use of effective stress and strain 


along with the potential function allows for the generation of master curves 
for elastic/plastic and elastic/viscoplastic material constants. 

If we let 

7 = ^<7 (13) 

then use of the overstress function ($) allows equation (6) to be written as: 

( 14 ) 


<*(»)> 


In terms of a functional relationship, we can use a power law expression 
to model test data and write the effective plastic strain rate as: 


F =< $(//) >= ( 


H 

IK 


(l/m) 


(15) 


where K and rn arc material constants found from test data on off-axis 
uniaxial specimens. 

Further examination of equation (14) leads to the following: 


7 < *(//) >= Sfju^u 

In terms of the in- plane components, this can be expanded to give: 

0 0 


(16) 


3 

2*< 


0 

<*n 

2a66 cr 12 J 


< 4>(//) > - 


9 (W 

\a 2 da 


0 <7^ 


0 

0 


0 0 \ a? u -ah 



i 

< 



\ *12 J 


(17) 


This allows the compliance matrix for the multiaxial plastic strain rate con- 
stitutive relation to be written as: 


3 




//■ 


c 

o 

0 ' 


0 T 22 * 

0 

.K. 

i 

o ■ 

°s 

?<ICG<?12 

^12 


(IS) 


Working at the lamina level and assuming plane stress conditions for an 
off-axis type of test under a state of uniaxial tension or compression (< 7 X ), 
the values of the principal stress components are: 


<7,, = <7 X cos 2 fl, <7-22 = a x sill* 0, (T\? = — <7 X SHI 0 COS 0 (19) 
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where 0 is the angle of the fiber direction relative to the load direction. Sun 
and Chen [3] showed that a function which depends on 0 and egg can be 
formed and used to help define effective stress and strain quantities. This 
function was defined as: 


"3 ~ 

h(0) = ^ (sin 4 0 + 2a m sin I 2 0 cos 2 

Using h(0), the effective plastic strain rate can be given as: 


t — 


fV 


m 


( 20 ) 


( 21 ) 


The overstress can also be written in terms of the axial stress components. 


II = h(0)(<7, - a‘ x ) (22) 

where < 7 * is the quasistatic stress and a x is the instantaneous applied stress. 
The quasistatic stress a* can be found by solving Sun and Chen’s [3] rate 
independent quasistatic constitutive relation: 


= £ + h(0)l n+ ')A(<r*) n (23) 

* J X 

where A and n are material constants found from test data. 

We note that the conditional statement on the overstress function ($) 
can now be written as: 


F =< $(//) >= 


$ if <7 r > <?* 
0 if a x < cr* 


(24) 


Expanding 4> we have: 




- <y*) 
K 


(l/m) 


(25) 


I his implies the axial plastic st rain rate can now ho completely written. 





(26) 
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If we let the quantity (3 be given by; 


? 


i / I \ l / m 

/)=[/.(«)]<'+*>(-) (27) 

we have a compact form for the axial plastic strain rate. 

e; = «<r, - <)'/" (28) 


For the time dependent inelastic behavior, we begin by utilizing the gen- 
eral multiaxial form; 

c' n = F(a,c,signd)a (29) 

noting that this implies the following relationships. 


F((T, e, signer) 


F(a,t) if a > 0 
0 if ct < 0 or <7 < <7* 


(30) 


Assume that F(a } e) comes from the quasistatic elastic/plastic expressions. 
We can make use of the quasistatic plastic compliance matrix [3]: 


[SP = 


0 0 
0 a\ 2 
0 0 


0 

0 


where the term 'P can be written: 

9 dl qp 
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$ = — ^ — = -nAa {n ~ 3) 
4 a 2 da 4 

Performing the necessary differentiation gives; 

d( qp 

<i ( , | 


and 


dt 

dc 22 da 


= 0 


+ 


(k qp 2 da ■ 


22 


dj%_ 

dl da dt ’ Da 2 2 dt 

^ - 3)^ + 

(it 'j(T 


(31) 


(32) 

(33) 

(34) 

(35) 
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and 


( 36 ) 


dt 

T 2 , 


*15 5cj2 da 0 c g t 2 da 12 


da dt da J 2 d< 

,^2 


dt = <7g(« - 3)a«,# — + 2'Pa 66 (T 1 2d 1 2 


(37) 


Collecting terms gives the function F as; 


a(n — 3) 


where the effective stress rate is: 


' 0 

0 

0 ' 


' 0 

0 

0 

0 

<722 

0 

+ # 

0 

oil 

0 

0 

0 

<Tj2 


0 

0 

^ a 66 a }2 . 


(38) 


<7 ~ ~{^22^22 + 2a^(T\2^l 2 ) 


(39) 


It should be noted that the quasistatic plastic compliance matrix [6'] ,p ap- 
pears as the second term in equation (38). 

Using equation (38), the multiaxial constitutive relation can be written 
in a more compact form as; 


0 "} = 


<r(n — 3) 
3i 


[.*T {*} + PT {*} 


(40) 


As before, assuming the case of uniaxial loading of an ofT-axis laminate, 
the inelastic strain rate can be written as: 

t'x = F(<r x , c n signer x )er x (4 1 ) 

Using the form of the uniaxial quasistatic elastic/plastic expression, we have; 

4” = ;2{h(0))}l n + l 'An(a x ) {n -"a x (42) 


This allows the F to be given in a specific sense as; 
F(a x , e x ,sign& x ) = | 


( .2[/i(fl)]^ n+l M7t(<r r ) (n-1 l if a x > 0 (loading) 

0 if <r x < 0 (unloading) 

or a,. < a* 


(43) 

where the multiplying factor of .2 gave good correspondence to test data for 
the material systems investigated. 
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CONSTITUTIVE MODEL, UNIAXIAL FORM 


Considering uniaxial loading of an ofT-axis laminate, it was assumed that 
the rate dependent strain could be decomposed into elastic and viscoplastic 
components. 


= c e + r v 

x 1 r 


('ll) 


Where the elastic term is; 



(«) 


and the viscoplastic term can be written in a general sense as; 


r'V _ 


f.r 

X 


+ e 


rP 


(46) 


Using the above expressions, four cases can be given to account for dif- 
ferent aspects of a test regime. Specifically, for tension loading, these cases 
are; 


1. Quasistatic; 


given; a x = < 7 * => h(0) < $ >= 0 and & x = 0 

we have; e x p = 0, i™ = 0, e x = ~ + h(0) n+l A(a x ) n 

* J X 

2. Loading; 


given; i x > 0 

we have; a x 


=► i x = h{6 ) < <P > +<j x + F 

~ \(i + -2[W] n+ Mu(a x )n-i)J 


3. Stress Relaxation; 


given; f x = 0 


<r,. < 0 


F = 0 


we have; = /t(0) <<!>>= or; h(0) < $ >= 


E x 
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4. Creep; 


given; a x = 0, i x > 0 => F = 0, e* = 0 

we have; i x =< f3(cr x — cr*p m > 

It is noted that the expressions for creep will not allow for prediction of creep 
recovery after unloading below the quasistatic stress level. 


LAMINATION THEORY 


The model presented above can be used to describe the nonlinear, rate de- 
pendent behavior in a laminated composite. For such a laminated composite, 
we can use typical notation from lamination theory so that Q<j is the stiff- 
ness matrix and S tJ is the corresponding compliance matrix. We also note 
that the transformed stiffness and compliance matrices are given by Qij and 
Sij respectively. This gives the constitutive equations for the quasistatic 
elastic/plastic case to be; 


{da}k = [Q]l{de e }k elastic (47) 

{<la}k = [Q]l P {dc'' p }k plastic (48) 

The subscript k refers to the individual or k ’ t, h layer in the laminate. If we 
let 

[OP = [QY + [OP (49) 

then, the combined expressions may be written for the quasistatic elas- 
tic/plastic case as; 


f (lN ] 


a i ir 

rqp 

<h° 

1 

\ 

► =r 


< 


{ dM J 


D | 1) 


{ (Ik _ 


(50) 


Similarly, for the elastic/viscoplastic case, we can write; 


[Or = [or + ion 


(51) 



then, the combined expressions may be written for the elastic/viscoplastic 
case as; 


( dN 1 
dt 


' A 

n ‘ 

CV]> 

f¥ 1 

- 

► = 

- 

~ 

i 


dM 
dt > 


.D | 

D 


i £ 1 


MATERIALS TESTING 

Material constants and related functions needed by this model were found by 
performing uniaxial tension or compression tests on off-axis laminates under 
isothermal conditions. By using effective values and collapsing this data into 
master curves, five experimentally derived constants were incorporated into 
the constitutive model. Additional information on the test procedures and 
equipment can be found in reference [9]. 

Test Specimens and Material 

Two polymer matrix composite material systems were* investigated in this 
study. The first, a graphite/thermoplastic was composed of Hercules 1 IM7 
fiber and Amoco 1 8320 matrix. The second material under study was a 
graphite/bismaleimide composed of Hercules IM7 fibers and Narmco 1 5260 
matrix. Both material systems had glass transition temperatures ( T g ) listed 
by the manufacture to be approximately 220°C'. 

Rectangular test specimens similar to those described in ASTM specifi- 
cation D3039-76 were cut from the finished panels. For the elastic material 
constants, tests were run on [0°]i2, [90°] 12 and [±45°]2 S specimens in order to 
determine E\, i/ l2 , E 2 and G x2 . Procedures similar to those outlined in ASTM 
specifications D3039-76 and D3518-76 were used to generate the constants. 

For the elastic/ viscoplastic material constants, off-axis tests on 15°, 30° 
and 40° coupons were performed using the rectangular specimen geometry 
described above. The specimens had an aspect ratio of 9.5:1 measured be- 
tween the specimen ends. Figure 1 shows a photograph of the test apparatus. 

All tests were conducted under isothermal conditions using either mono- 
tonic tension or compression loads. The four temperatures selected for study 

'The use of trade names in this paper docs not constitute endorsement, either expressed 
or implied, by the National Aeronautics and Space Administration. 
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were 23°, 70°, 125° and 200°C. The tests were run using strain or load control 
as appropriate. 

Material Constants 

Five material constants are required by the analytical model for any given 
temperature and load direction. These constants are: a m for the potential 
function, A and n for the quasistatic elastic/plastic relations, and K and m 
for the rate dependent elastic/viscoplastic relations. These constants were 
all found using data from simple off-axis tension or compression tests. 

As shown in reference [9], using a strain controlled test and an approach 
similar to that outlined by Yen [10], all of the constants were extracted from 
uniaxial tests with repeated holds built into the test to allow for stress re- 
laxation to occur. The idea behind this approach is that during the period 
of stress relaxation, the stress will decrease rapidly towards some limiting 
value. This limiting value is assumed to be the quasistatic stress and repre- 
sents the stress needed to solve the elastic/plastic constitutive equation for 
a given strain. By repeating these periods of stress relaxation during the 
course of the test, enough quasistatic points are found to allow an effective 
stress, effective plastic strain curve to be constructed. This curve represents 
the quasistatic behavior. The effective curves from different off-axis tests can 
then be plotted together and collapsed into a single curve by selecting the 
appropriate value of a f>6 . A power law fit to this curve determines A and n. 

During stress relaxation, the total strain rate is zero, therefore, as shown 
by the expressions in case 3, the viscoplastic term can be equated to the 
stress rate divided by the elastic modulus. During relaxation, the overstress 
and the stress rate can be defined for any point in time. This allows the 
construction of a overstress (//) versus effective viscoplastic strain rate ($) 
graph. These master curves are fit with a power law expression which gives 
the constants K and m. 



RESULTS 


Not all specimens were tested to failure, so failure levels are not reported. 
Optical observations and edge replication revealed no matrix cracking prior 
to specimen failure. 

A review of material properties is provided in table 1, which also gives a 
complete list ot all the constants. The material constant a gg which appears in 
the potential function was found to be independent of temperature. However, 
it was found to be dependent upon material system. In addition, a m for 
IM7/8320 was found to be dependent on whether the material was in tension 
or compression. 

From the quasistatic master curves such as those shown in figures 2-3, it is 
apparent that the IM7/5260 material shows less tendency towards nonlinear 
elastic/plastic behavior than the IM7/8320 system. Both systems show a 
definite trend towards increased ductility as temperature increases. 

From the rate dependent master curves such as shown in figures 4-5, one 
can see from a comparison of the two systems that the IM7/5260 material 
shows a greater tendency towards higher overstress than the IM7/8320 does 
for an equivalent increase in plastic strain rate. This trend implies that the 
IM7/5260 is exhibiting less viscoplastic behavior than the IM7/8320. Both 
material systems indicate a slight increase in the viscoplastic behavior with 
an increase in temperature. 

For both material systems, the master curves for tension dilTer than those 
for compression at the same temperature. This may imply that the defor- 
mation mechanisms in tension and compression are not the same. 

On the lamina level, the. uniaxial elastic/viscoplastic constitutive model 
was solved numerically using the fourth order Rungc-Kutta technique with 
a modified Newton technique to find the roots of the quasistatic equation. 
Comparing the uniaxial test results against the analytical model shows that 
the model does well in predicting several phenomenon including: stress re- 
laxation, sliorl term creep, linear elastic unloading and variable strain rate 
loading. Typi< al comparisons between test and predictions are given in fig- 
ures 6-9. 

Short term creep behavior of a 15° olT-axis specimen at 70°C is shown in 
figure 6. The applied stress history, resultant strain history and predicted 
strain history are given. The correlation between test and predicted strain 
is typical for short term creep tests performed on both material systems. 
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The prediction of creep behavior is a good verification of the model since the 
material constants used for the creep prediction were found from the stress 
relaxation procedures described previously. 

Figure 7 shows the relationships between stress, strain and time for a 
25° olf-axis tension specimen at 70°C. In this case, the test was run under 
strain control and the resultant stress was measured. Periods of constant 
strain rate, stress relaxation and unloading comprise the input strain history. 
The predicted stress/time and stress/strain behavior is plotted against the 
measured values in figures 7a and 7b. A good correlation between test and 
predicted values is evident. Figure 8 is similar to figure 7 except for the 
higher test temperature. Once again, good correlation is shown. It should 
be noted that the 25° data was not used to construct the master curves and 
therefore correlation between test and prediction can be used to help verify 
the model. 

Figure 9 gives test and predicted values for a strain controlled, uniaxial 
compression test of a 30° olf-axis specimen at 125°C. The results indicate 
a reasonable correlation, but as was typical with most of the compression 
results, the comparison between test and predicted values was not as good 
as the tension cases. 

In general, comparison of test versus predicted values gave good results 
for both materials, in tension and compression, and throughout the range of 
test temperatures. 


SUMMARY 

An elastic/viscoplastic constitutive model was developed to describe the ob- 
served nonlinear, rate-dependent behavior in advanced polymer matrix com- 
posites at elevated temperatures. Formulations for the general multiaxial 
case and the specific uniaxial case were given. Both uniaxial tension and 
compression loading of off-axis specimens were investigated experimentally. 

Based upon the need for material constants in tin* clastic/viscoplastic 
constitutive model outlined above, test procedures were developed for the 
determination of the required parameters, 'flic test methods used were found 
to be reliable and repeatable. Simple ofT-axis tests with repeated periods 
of stress relaxation were found to be useful for generating the necessary 
data. The specimen geometry and end tabs minimized the extension-shear 
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coupling, however, test temperatures close to T g will enhance this behavior. 

Comparisons between test data and predicted values indicated that the 
model provides a reasonable prediction of the behavior of load or strain con- 
trolled tests. Periods of loading, unloading, stress relaxation and creep were 
accounted for. Observed material behavior that the model does not account 
for include: creep recovery after unloading and nonlinear time dependent 
behavior during unloading. A complete description of cyclic behavior would 
require that these phenomenon be correctly modeled. 

It is expected that further work with this constitutive model will sup- 
port the development of a laminate level computer code which can be used 
to predict the elastic/viscoplastic behavior of any arbitrary laminate under 
conditions of combined mechanical and thermal loads. Further testing on 
IM7/5260 and IM7/8320 laminates should provide data for additional verifi- 
cation of the analysis procedures. 
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Material 

Type 

°C 

Elastic 

E 1 (GPa) | E 2 (GPa) | G u (GPa) 

/y 12 

Elastic/Plastic 
age | A(MPa)~ n | n 

Elastic/ Viscoplastic 
K(MPa) I m 

IM 7/5260 
Tension 

23 

152.8 

8.7 

5.2 

.30 

.60 

5.07E-10 

IHa 

lUlilaUlkl 

.92 

70 

161.7 

9.2 

5.8 

.31 

.60 


iiftfci 

2.09E+05 

.98 

125 


8.8 

5.2 

.36 

.60 

1.18E-09 


HMD'H 

1.03 

200 

■HIM 


5.1 

.35 

.60 


EMi 

9.83E+04 

.94 

IM7/5260 

Compress. 

23 

152.8 

8.7 

5.2 * 

.30 

.60 

1.59E-09 

3.02 

1.24E+06 

1.12 

70 

161.7 

9.2 

5.8 * 

.31 

.60 




.83 

125 

156.5 

8.8 

5.2 * 

.36 

.60 



2.59E+05 

1.01 

200 

154.3 

7.5 

Cn 

* 

.35 

.60 

4.30E-14 

6.91 

9.73E+04 

.92 

IM7/8320 

Tension 

23 

157.9 

7.1 

4.3 

sai 

.30 




.79 

70 

153.8 

7.9 

4.3 

.34 

m 

2.19E-08 

3.36 

I3.33E+04 

.86 

125 

142.0 

7.5 

3.2 

.35 

mm 

3.61 E- 11 



.72 

200 
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.83 
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.35 | 

m 
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Table 1: Material Properties and Constants. * Assumed from tension data 
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Figure 2: Elastic/plastic master curve's for tension loading. 
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Figure 3: Elastic/plastic master curves for compression loading. 
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Figure 5: Elastic/viscoplastic master curves for compression loading. 
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Figure 6: Comparison of test results and model predictions for creep in 
tension, of a typical IM7/8320 15° off-axis specimen at 70°C under load 
control. 
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Figure 7: Comparison of test results and model predictions of a. typical 
IM7/5260 25° off-axis tension specimen at 70° C under strain control. 
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Figure 8: Comparison of test results and model predictions of a typical 
IM7/5260 2.'3° off-axis tension specimen at 125°C under strain control. 
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Figure 9: Comparison of test results and model predictions of a typical 
1M7/5260 30° ofF-axis compression specimen at 125°C under strain control. 
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